Acute and Chronic Changes in Gene Expression After CMV DNAemia in Kidney Transplant Recipients

Cytomegalovirus (CMV) viremia continues to cause significant morbidity and mortality in kidney transplant patients with clinical complications including organ rejection and death. Whole blood gene expression dynamics in CMV viremic patients from onset of DNAemia through convalescence has not been well studied to date in humans. To evaluate how CMV infection impacts whole blood leukocyte gene expression over time, we evaluated a matched cohort of 62 kidney transplant recipients with and without CMV DNAemia using blood samples collected at multiple time points during the 12-month period after transplant. While transcriptomic differences were minimal at baseline between DNAemic and non-DNAemic patients, hundreds of genes were differentially expressed at the long-term timepoint, including genes enriching for pathways important for macrophages, interferon, and IL-8 signaling. Amongst patients with CMV DNAemia, the greatest amount of transcriptomic change occurred between baseline and 1-week post-DNAemia, with increase in pathways for interferon signaling and cytotoxic T cell function. Time-course gene set analysis of these differentially expressed genes revealed that most of the enriched pathways had a significant time-trend. While many pathways that were significantly down- or upregulated at 1 week returned to baseline-like levels, we noted that several pathways important in adaptive and innate cell function remained upregulated at the long-term timepoint after resolution of CMV DNAemia. Differential expression analysis and time-course gene set analysis revealed the dynamics of genes and pathways involved in the immune response to CMV DNAemia in kidney transplant patients. Understanding transcriptional changes caused by CMV DNAemia may identify the mechanism behind patient vulnerability to CMV reactivation and increased risk of rejection in transplant recipients and suggest protective strategies to counter the negative immunologic impact of CMV. These findings provide a framework to identify immune correlates for risk assessment and guiding need for extending antiviral prophylaxis.


INTRODUCTION
Despite years of improvements in screening, prophylaxis, and treatment options, CMV continues to cause significant negative impact in solid organ transplant recipients, via both direct and indirect effects (1). Direct effects include disseminated infection and end-organ disease including pneumonitis, esophagitis, and colitis. Indirect effects, in contrast, cause more impactful consequences of CMV infection, including vulnerability to secondary infections and development of alloimmunity, leading to allograft dysfunction. Heterologous immunity has been postulated as the mechanism behind this negative effect on the allograft, although generalized increase in inflammation is another possible cause (2,3). Evaluation of the transcriptional changes secondary to CMV primary infection or reactivation of CMV can provide insight into the mechanisms behind allograft dysfunction secondary to CMV.
Previous studies on host response to CMV are primarily derived from the MCMV mouse model. These have demonstrated upregulation of Cxc3r1, Cd69, and Cd103, associated with inflation of T cell memory, tissue association, and exhaustion, as well as transcription factors associated with inflammation (4,5). This shift towards CMV-specific memory T cells is postulated to underlie the progression to increased immune senescence (6). Another key aspect of CMV infection is the establishment of latent infection and virus-mediated downregulation of host genes, blocking ability to clear the virus (7). However, less is known in the context of human CMV infection.
Solid organ transplantation presents a unique clinical model to investigate CMV infection in humans, where the serostatus of both donor and recipient is known prior to transplantation, and CMV PCR is followed for surveillance purposes after transplantation, allowing for identification of the initiation of primary CMV infection or reactivation. This provides an opportunity to examine changes in host gene expression before and after the start of detectable CMV DNAemia, providing insights into the impact of CMV infection on the transplant recipient's innate and adaptive antiviral immune response and potential influence on graft outcome. We therefore characterized the whole blood transcriptome from a cohort of kidney transplant recipients over several time points during and after detectable CMV infection to analyze changes in gene expression over time, as well as differences between patients with or without CMV DNAemia.

Patient Recruitment and Clinical Care
Kidney transplant recipients were enrolled in a UCLA IRBapproved study, and blood was collected at regular intervals after transplantation followed by PBMC isolation and storage, as previously described (8). Patients received induction with either anti-thymocyte globulin (ATG) or basiliximab depending on pre-transplant levels of sensitization and donor kidney quality followed by protocol-based immunosuppression with tacrolimus, mycophenolate mofetil, and prednisone, as previously described (8). Patients intolerant of tacrolimus were switched to cyclosporine. CMV prevention was performed according to local protocols as previously described, summarized as follows: 6 months of valganciclovir for highrisk donor positive (D+) and recipient negative (R−) patients and 3 months of valganciclovir for intermediate-risk recipient positive (R+) patients who received ATG induction. R+ patients who received basiliximab or low-risk (D−/R−) patients received acyclovir prophylaxis to prevent HSV and VZV infection. All patients regardless of prophylaxis regimen underwent regular CMV PCR screening to detect CMV DNA in peripheral blood for the first year after transplantation.

Patient Selection
Patients with a history of positive testing by CMV DNA by PCR were identified by review of microbiology records. We reviewed our repository database to determine whether PBMCs were available that corresponded to the following timepoints: (1) Baseline (prior to DNAemia start), (2) 1 week post-DNAemia start (week 1), (3) 1 month post-DNAemia (month 1), and (4)~1 year after transplantation (long-term). Research blood was collected from stable outpatients at the time of kidney transplant clinic visit. Four of the 31 patients that developed CMV DNAemia after transplantation were donor IgG seropositive, recipient seronegative; the others were recipient IgG seropositive. For patients with multiple episodes of CMV DNAemia over 137 IU/ml, the first episode closest to transplantation was studied.
Patients with history of CMV DNAemia were matched on a 1:1 basis to a cohort of kidney transplant recipients without history of CMV DNAemia based on deceased versus living donor status, patient age, sex, race/ethnicity, and induction type with kidney transplant recipients who were either D+/R− or R+. Samples were selected for each control patient that corresponded in terms of time post-transplant with the baseline and long-term CMV DNAemia patient samples. Of these patients negative for CMV DNAemia, 24 were CMV seropositive and 7 were seronegative with CMV-positive donors. Antiviral prophylaxis was administered by protocol, so that high-risk patients (D+/R−) received 6 months of Valcyte, and intermediate-risk patients (R+) received 3 months of Valcyte if they received antithymocyte globulin induction, otherwise acyclovir for 3 months, as previously described (9). Patients with detectable levels of CMV DNAemia >137 received antiviral therapy per protocol for a minimum of 2 weeks, and until CMV PCR was <137, per American Society of Transplantation guidelines (1). We focused exclusively on CMV infection as this was our central scientific question, and due to the fact that other post-transplant infections are either much rarer compared with CMV, the most common viral infection after kidney transplantation (10), or site specific and limited to a single organ such as the case for urinary tract infection, and therefore predicted to have a much less significant impact on peripheral blood transcriptional analysis compared with CMV.

Sample Preparation and Sequencing
RNA was isolated from 1 ml RBC-lysed whole blood samples stored in RNAlater (500 ml). Globin RNA was removed using the GLOBINCLEAR kit. Libraries were prepared for samples that passed quality control using the KAPA stranded mRNA kit. Single-read sequencing with read length of 50 bp was performed using the Illumina HiSeq3000 platform.

Read Quality Control and Alignment
After sequencing, quality control of reads was checked using FastQC 0.11.8 (11). Reads were aligned to the GRCh38 human reference genome using STAR 2.4.2 (12). We used samtools (13) and picard (https://github.com/broadinstitute/picard) to index the resulting BAM files and remove any PCR duplicates. Ten samples were removed from further downstream analysis because of low library size and/or a high proportion of reads that did not uniquely map to the reference. We used featureCounts (14) to produce the final count table.

Univariate Statistical Analysis
To test for association between demographic and clinical data and CMV DNAemia status, Fisher's exact test was implemented for categorical variables, while Welch's two-sample t-test was implemented for continuous variables.

Differential Expression Analysis and Pathway Analysis
We used the edgeR 3.9 package (15) in R (16) to convert the raw counts to counts-per-million (CPM), to perform trimmed mean of M-values normalization, and to test for differential expression by fitting a negative binomial generalized linear model, estimate dispersion estimates, and perform a likelihood ratio test. We tested for differential expression between (1) baseline and week 1, (2) baseline and month 1, (3) baseline and long-term, (4) CMV DNAemia vs no CMV DNAemia at baseline, and (5) CMV DNAemia vs no CMV DNAemia at long-term. We used Ingenuity Pathway Analysis (IPA) to perform canonical pathway enrichment analysis of the differentially expressed genes (DEGs).

Time-Course Gene Set Analysis
We implemented the TcGSA R package (17) to perform timecourse gene set analysis to longitudinally test the stability of gene set expression. In summary, for each gene set, TcGSA implements a linear mixed model of gene expression that includes a fixed-effect for the expression of each gene in a gene set, a random-effect of the expression of each gene in a gene set from each patient, as well as a time function to model the time trend for each gene. This time function can be divided into a fixed effect that represents the average trend of the gene set and a random effect that accounts for the heterogeneity that may be present between each gene and the gene set and may take a linear polynomial, cubic polynomial, or natural cubic spline form. For this study, a cubic polynomial form was chosen for the time function as it provided a better fit to the data than a linear polynomial form but an equivalent fit to the natural cubic spline form but with fewer degrees of freedom. After modeling the gene expression in each gene set, a likelihood ratio test was implemented to test whether a gene set had a significant time trend, where the null hypothesis was that the expression genes in a gene set were stable over time.

Patient Demographics
A summary of patient demographics is provided in Table 1. As expected, due the propensity-matching design for patient selection, age, sex, race, and induction type were comparable between patients with DNAemia (≥137 IU/ml of CMV detected in a PCR test; abbreviated PCR+) and non-DNAemic patients (PCR−) ( Table 1). A similar number of patients were high risk (D+/R−) as intermediate risk (R+) in both groups, and median GFR at 3 months and rates of rejection were also similar between groups. Patients experiencing DNAemia were PCR-positive at a median of 80 days after transplantation (range: 10 to 561 days). Median peak viral load was 757 IU/ml (range 146 to 13,900 IU/ ml). Two patients with DNAemia were diagnosed with CMV syndrome or CMV end-organ disease based on standard definitions (18).

Transcriptomics of DNAemic Patients Over Time: Individual Gene Analysis
To characterize the transcriptional landscape of CMV DNAemia patients across a 1-year timeframe post-transplant, we performed RNA-seq on whole blood from 31 patients with CMV DNAemia, across the four timepoints (n=131) described above. Both differential expression analysis and pathway analysis were performed for each time point.
Review of differential expression of all genes with FDR<0.1 demonstrated striking changes after detection of CMV DNAemia. Analysis of the difference between baseline and week 1 post-DNAemia revealed the most prominent transcriptional changes with respect to the number of DEGs, with 2,456 DEGs total and 538 of those DEGs having an absolute log fold-change (|logFC|) ≥ 1 at FDR ≤ 0.1 ( Figure 1A; Table 2A; see Supplementary Data 1 for full listing of all DEGs). A total of 1,299 DEGs were upregulated and 1,157 downregulated at week 1 after CMV DNAemia compared with baseline ( Table 2A). Several of these DEGs encompass genes that are involved in immune signaling, cytotoxicity, and T cell activation, such as TRGV2 (T Cell Receptor Gamma Variable 2) and TRDV2 (T Cell Receptor Delta Variable 2), which are upregulated, and IL1A (Interleukin 1 Alpha), which is downregulated (top 20 up-/downregulated genes are shown in Table 2B). Besides these notably higher transcriptional changes, we detected subtle variations in several immunologically relevant genes at week 1 post-viremia, including transcripts encoding for transcription factors HOPX (Homeodomain-only protein) and BATF2 (Basic Leucine Zipper ATF-Like Transcription Factor 2) were slightly upregulated 1 week after CMV DNAemia (Supplementary Data 1). In support of potential ongoing TCR engagement, we found an increase of CD69 and a simultaneous expression of TCR dependent co-inhibitory molecules of LAG3 (Lymphocyte Activating 3) and CTLA4 (Cytotoxic T-lymphocyte Associated Protein 4) at 1 week post-DNAemia. In addition, we found that transcripts encoding DUSP4 (Dual Specificity Phosphatase 4), KLRG1 (Killer Cell Lectin Like Receptor G1), CD70 were also enriched, which are linked to inhibitory function and to negative regulation of T cell activation and proliferation.
While many of these genes a showed transient change in expression, several molecules showed a stable pattern of expression emerging week 1 and remained detectably higher at

A B C
FIGURE 1 | Volcano plots of genes that are differentially expressed between (A) baseline and week 1; (B) baseline and month 1; (C) baseline and long-term in patients that developed DNAemia. Gray points represent genes that are not differentially expressed, green points represent genes that have a log 2 FC ≥ 2 but not differentially expressed, blue points represent genes that are differentially expressed (FDR ≤ 0.1) but have a log 2 FC < 2, and red points represent genes that are both differentially expressed (FDR ≤ 0.1) and have a log 2 FC ≥ 2.
the month 1 and long-term timepoints. One such molecule is FCRL6 (Fc Receptor Like 6) (Table 2B), a distinct marker of cytotoxic effector lymphocytes. The other gene is GPR56, a marker associated with CD56null/CD16+ NK cells. Interestingly, we also noted decrease in IRF8, which orchestrates adaptive NK-cell responsiveness and antiviral immunity.
Analysis of differentially expressed genes at 1 month after CMV DNAemia compared with baseline revealed ongoing changes in pathways important for immunologic function, with a total of 1,972 significantly enriched DEGs ( Figure 1B; Table 2A). These included upregulation of genes important to adaptive and innate immune function including TRGV2, TRDV3   (T Cell Receptor Delta Variable 3), IGFL2 (IGF Like Family Member 2), CCL4 (C-C Motif Chemokine Ligand 4), and IFNG (Interferon-g) ( Table 2B).
At the long-term timepoint, there were fewer changes in canonical pathways compared with the week 1 to baseline change; nonetheless, 29 statistically significant differences in DEGs were detected ( Figure 1C; Table 2A). Genes from within these DEGs included upregulation in TRGV2, GZMH, GZMA, FCRL6 (Fc Receptor Like 6), and CCL5 (Chemokine Ligand 5) ( Table 2B). Although it was not possible to clearly differentiate the direction of change and assign z-score for the long-term timepoint, several important pathways were identified as undergoing significant change after resolution of DNAemia ( Figure 2; see Supplementary Datas 2 and 3 for a full listing of pathways and genes for each pathway). This timepoint was notable for changes in pathways important for innate immune and cytotoxic T cell function, including Communication between Innate and Adaptive Immune cells (-log(p) = 2.39), Granulocyte Adhesion and Diapedesis (-log(p) = 1.86), and Granzyme A Signaling (-log(p) = 1.7) (Figure 2). Long-term impact on T cell function is reflected by the several pathways related to IL-17 signaling, as well as relative changes in the Hematopoiesis from Pluripotent Stem Cells (-log(p) = 1.32) and Primary Immunodeficiency Signaling (-log(p) = 1.31) pathways. These results may suggest that in the long-term after resolution of CMV DNAemia, both the innate and adaptive immune responses are impacted.

Evaluation of Differentially Expressed Genes at 1 Week Post-DNAemia Across All Timepoints
To further visualize longitudinal changes in gene expression post-CMV DNAemia, we generated a heatmap of the logCPM expression of the DEGs between baseline and week 1 in PCR+ patients with |logFC| ≥ 2 mapped across the four timepoints (baseline, week 1, month 1, and long-term) (Figure 3). We observed multiple visible trends of expression, with many genes becoming either highly up-or downregulated at week 1. There is a striking difference between baseline and week 1, with multiple genes upregulated at week 1 including FCGR1CP, NCPAH, and KIF2C. Interestingly, at month 1 post-viremia, this upregulation is beginning to decrease, so that there is a return to near-baseline levels at the long-term sample although notably, as with the pathway analysis, expression levels for

Time-Course Gene Set Analysis
To determine which genes demonstrated persistent longitudinal changes, we evaluated time-trends across all four time points using TcGSA. We found that 15 of the 18 canonical pathways from baseline versus week 1 had a statistically significant time-trend across all four timepoints (p ≤ 0.05) ( Figure 4A; Table 3). We observed two major patterns of gene expression: (1) downregulation of a pathway at the week 1 timepoint, followed by a gradual upregulation to the long-term timepoint or (2) upregulation at week 1 followed by a gradual downregulation to the long-term timepoint. We also observed that most of the pathways contain genes that exhibit both patterns, which leads to there being more than one significant time-trend detected by TcGSA. For example, CD27 Signaling in Lymphocytes and Oxidative Phosphorylation show two timetrends, while Interferon Signaling only has one time-trend ( Figure 4B). We did not observe any significant time-trends for any of the canonical pathways identified as changing significantly from week 1 versus month 1 or month 1 versus long-term timepoints.

Transcriptional Differences Between At-Risk Patients With and Without CMV DNAemia
To compare the transcriptional profile of patients who developed CMV DNAemia with those that did not develop DNAemia posttransplant, we also performed RNA-seq on 31 matched kidney transplant recipients without CMV DNAemia. We carried out differential expression tests between the following sets: (1) DNAemia versus non-DNAemia at baseline and (2) DNAemia versus non-DNAemia at the long-term timepoint.
We identified a total of 16 DEGs at FDR ≤ 0.1 with six of those DEGs having an |logFC| ≥ 1 between patients who did or did not develop DNAemia at the baseline timepoint, prior to clinical detection of CMV by PCR testing ( Table 4; see Supplementary Data 1 for a full differential expression analysis results). One of the top upregulated genes for CMV DNAemia patients at baseline is USP18, a deubiquitinating protease that may play a role in interferon responsiveness. Interestingly, we could also detect PDE5A, associated with modulation of chronic inflammation derived myeloid-derived suppressor cells and regulatory T cells activity. We also noticed significant upregulation of PPP1R10 (protein phosphatase 1 regulatory subunit 10) and the transcription factor NFAT5, which regulated type 1 interferon responses, within the CMV DNAemia group at baseline. However, SIGIRR (Singleimmunoglobulin interleukin-1 receptor-related), an inhibitor of TLR and IL-1R signaling, showed significant downregulation.
In contrast to the results from the baseline timepoint, there were a greater number of statistically significant differences observed at the long-term timepoint. There were 1,420 DEGs between patients with and without history of CMV DNAemia at the long-term timepoint at FDR ≤ 0.1 ( Table 4 and Figure 5A). We observed 21 canonical pathways that were significantly enriched for (p ≤ 0.05, |z-score| ≥ 2) including multiple pathways important in innate cell functioning. These included Fcg Receptor-mediated Phagocytosis in Macrophages and Monocytes (-log(p) = 7.84, z-score = 4.02), Interferon Signaling (-log(p) = 6.22, z-score = 3.67), and IL-8 Signaling (-log(p) = 3.71, z-score = 3.67) ( Figure 5B; see Supplementary Data 2 and C for a full listing of pathways and genes for each pathway). Interestingly, all immunologic IPA pathways demonstrating significant changes were in the direction of upregulation in patients with history of DNAemia compared with those without, which persisted despite the fact that for all patients in this cohort, DNAemia had resolved at the time of the long-term timepoint.

DISCUSSION
This study surveyed the transcriptomic landscape of longitudinally collected peripheral blood leukocytes from kidney transplant recipients with and without CMV DNAemia.  Kidney transplantation represents a unique opportunity to study the time course of CMV reactivation or primary infection given the administration of immunosuppression and close clinical monitoring for CMV by PCR in peripheral blood.
A major finding of this study is the early upregulation of gene expression followed by long-term changes in gene transcription from peripheral blood cells in patients experiencing DNAemia post-transplantation, which are not observed in patients that do not experience CMV DNAemia. These changes occurred early after DNAemia, with significant impact on transcription as early as 1 week after DNAemia start, but also persisted following resolution of DNAemia. There were detectable differences in gene expression in patients with DNAemia at the long-term compared with the baseline time point (Figure 1). In addition, we detected baseline differences in gene expression in patients with history of DNAemia compared with those without history of DNAemia ( Figure 5). These changes include genes important in the innate immune response as well as T cell activation and cytotoxicity, confirmed at the individual gene level as well as via canonical pathway analysis.
The changes in gene expression were most striking at the 1week post-viremia timepoint as compared to the baseline sample. These included T cell specific transcripts TRGV2 and TCRDV2 important in T cell receptor expression including in gamma delta T cells, and LAG3, CTLA4, and KLRG1, associated with T cell exhaustion. At the 1-month time point, other genes important in adaptive and innate immune function were noted including CCL4 and IFNG, as well as FCLR6 and GPR56, which are associated with immune aging.
While it is not surprising that the biggest impact was seen at this early timepoint, several genes were persistently altered following CMV DNAemia, leading to development of a new baseline for gene expression after resolution of CMV DNAemia. Genes in pathways such as CCR5 Signaling in Macrophages and Cytotoxic T Lymphocyte-mediated Apoptosis of Target Cells remain upregulated even after resolution of DNAemia, as shown in Figure 2. Of note, Interferon Signaling and the STAT3 Pathway are pathways that are downregulated in patients experiencing CMV DNAemia, yet still upregulated 9-18 months post-DNAemia relative to the patients who never experienced CMV DNAemia. Many of the individual genes seen early post-CMV DNAemia, such as TRGV2, FCLR6, and CCL5, remain elevated at the long-term timepoint. In addition, PPP1R10 is known to play a role in cell cycle progression, DNA repair, and apoptosis by negatively regulating the activity of protein phosphatase 1, a proapoptotic activator important in DNA repair. This suggests that the experience of CMV DNAemia, even after resolution, leaves a lasting impact on the immune system, which promotes inflammation and immune dysregulation. Interestingly, while six patients did experience two or more episodes of CMV DNAemia, these subsequent events were generally shorter, had lower viral loads than the index event, and did not meaningfully affect gene expression profiles when DNAemia recurrence was adjusted for in differential expression analysis (data not shown).
A positive association between CMV infection or reactivation has previously been shown by epidemiologic studies demonstrating increased rates of rejection and graft loss in high-  risk as compared with intermediate-or low-risk kidney transplant recipients months or years after initial infection, despite receipt of antiviral prophylaxis (19). Moreover, CMV infection/reactivation can cause dysregulation in immune system (20,21). This imbalance in the immune system may increase the risk of transplant rejection (2). Notably, many of the genes identified in our analysis as upregulated at 1 month after CMV DNAemia begins have also been identified in association with antibodymediated or T cell-mediated rejection, including CCL4, IFNG, and TRDV3, as well as FGFBP2, which is elevated at the long-term timepoint compared with baseline timepoint ( Table 2B). Upregulation of these genes may increase inflammation and rejection risk in the transplant allograft, namely, chemokines and chemokine receptors, and cytokines such as IL-6 (22). Other gene pathways including those impacting CCR5 have been previously identified as important in pathogenesis of rejection and chemotaxis of immune cells to the allograft (23). This observation, therefore, suggests a model by which the CMV response leads to the generation of heterologous antigen-specific T cells with cross-reactive alloimmune specificity to donor HLA molecules that can mediate allograft injury (24)(25)(26). Thus, a polyclonal antigen-specific immune response to CMV induces pro-inflammatory transcriptional upregulation in peripheral blood cells that peak early after CMV infection, but even after resolution of DNAemia, there are long-term changes in gene expression in peripheral blood cells that cause a longstanding impact in the kidney allograft ( Figure 6A). Given that this association between CMV infection and rejection is seen in a variety of transplant recipient types (27), it is likely that this model may hold for other solid organ transplant recipients, and may explain the hypothesized connection between CMV and heart disease and other comorbidities (28)(29)(30)(31).
The other notable finding is the lasting changes in adaptive immune response, impacting T cell senescence and exhaustion likely promoted by the fact that after resolution of clinically detectable viremia, the subclinical or latent virus triggers persistent CMV antigen exposure ( Figure 6B). We hypothesize that in contrast to the pro-inflammatory effect of early DNAemia, this long-term effect is the likely mechanism behind vulnerability to other secondary infections including opportunistic bacterial and fungal infections that have been observed in transplant recipients (32)(33)(34). In our analysis, we noted upregulation of transcription factors associated with exhaustion including HOPX, which has been shown to be highly upregulated in terminally differentiated effector/memory Th1 cells and positively regulate effector differentiation, function, persistence, and survival of T cells (35). We additionally noted upregulation of Batf2, a transcription factor that belongs to the basic leucine zipper transcription factor family, which includes BATF and BATF3. BATF2 is expressed in immune cells such as T cells, B cells, macrophages, and dendritic cells and IFN-g-inducible BATF2 in innate immune cells controls Th17mediated immunopathology by suppressing IL-23 production, as well as FCRL6, which is upregulated in CD56dim NK cells, Vdelta1+ and Vdelta2+ gamma-delta T cells, effector and effector memory CD8+ T cells, and rare cytotoxic CD4+ T cells in diseases characterized by chronic immune stimulation (36). This pro-inflammatory process likely drives CMV and agingassociated T cell senescence, leading to premature increase in frequency of senescent and exhausted T cells in individuals with history of CMV exposure, even in healthy adults without endorgan disease or history of transplantation. Understanding the mechanisms behind the deleterious effects of CMV at different timepoints of infection can suggest possible interventions, and ideal timing of interventions, to block these negative outcomes. Many attempts have been made to establish risk stratification for CMV DNAemia as a tool in organ transplantation, including detection of a CMV-specific T cell response (37,38). We found that at the baseline sample prior to CMV DNAemia, there were already differences in gene expression between those who developed versus those who never develop CMV DNAemia ( Table 4; Supplementary Data 1). One of these genes, USP18, is related to interferon responsiveness (39), raising the question of whether this indicates increased vulnerability to CMV, as opposed to the detection of low levels of CMV replication below the limit of detection of the clinical PCR test. Further exploration of this observation may reveal more about the mechanism behind vulnerability to CMV, especially reactivation of CMV in CMV seropositive patients. In addition, it suggests the possibility of using monitoring of host transcription for early diagnosis and risk stratification of CMV infection or reactivation.
One limitation of this study is that most patients, in both the CMV DNAemia and no CMV DNAemia arms, have a history of CMV infection. Therefore, it is likely that the majority of the transcriptomic profiles of CMV DNAemic patients represent CMV reactivation and not primary infection. This issue also makes it difficult to compare the impact of high risk (D+/R−) status versus intermediate risk (R+) status on acute and longterm changes in gene expression. Another potential limitation of this study is the fact that all patients received antiviral therapy after CMV DNAemia was detected, while the patients without CMV DNAemia received antiviral prophylaxis and monitoring only. However, this limitation is mitigated by the fact that as a single-center study, all patients received similar regimens of immunosuppression, monitoring for rejection, antiviral prophylaxis, and CMV PCR monitoring. In addition, the most striking differences in gene expression were observed early after CMV DNAemia, likely before the start of antiviral therapies, which in addition would be unlikely to directly affect gene expression. We additionally lacked data on vaccination in this cohort, which would have largely consisted on influenza vaccine during this time period; however, we would predict that the impact of CMV viremia would overshadow any vaccination impact, which would have occurred in similar frequency in the CMV infection compared with the control patient cohorts. Another limitation of this study is the lack of associations between study participants and favorable or adverse clinical outcomes because of the relatively small size of our cohort and the low incidence of either favorable or adverse clinical events. We plan to perform future studies with a case-control study design that can properly address this limitation. Yet another limitation to this analysis is that as patients are followed longitudinally over the course of an episode of CMV DNAemia, they are also progressively farther away from induction immunosuppression, which may also have an impact. One meliorating factor, however, is that the control group is also followed over the year after transplantation, allowing analysis of changes related to CMV infection as opposed to change in time from transplantation alone.
Finally, as we performed bulk RNA-seq on whole PBMC samples, we are not able to characterize the longitudinal transcriptomic profiles of specific cell populations. In future studies, we plan to perform transcriptomic profiling of CMVtetramer positive T cells or single-cell RNA-seq to overcome this limitation.
These data provide new insights into in vivo dynamics of the longitudinal immune response to CMV DNAemia, which have the potential to identify mechanisms of susceptibility to CMV infection and allograft rejection.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi. nlm.nih.gov/geo/query/acc.cgi?acc=GSE168598.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by UCLA Institutional Review Board. The patients/ participants provided their written informed consent to participate in this study.

FUNDING
This study was supported by a National Institutes of Health grant U19 AI128913 (PIs: ER and MS). The funding bodies played no role in the collection, analysis, or interpretation of data in this study.